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SUMMARY 


The  purposes  of  this  research  program  were  to  extend  or  develop  analytical 
methods  for  determining  rotor  blade  aeroelastic  stability  limits  and  to 
perform  stability  calculations  over  a  range  of  design  and  operating  vari¬ 
ables  for  articulated  and  nonarticulated  configurations.  The  usefulness  of 
simpler  analytical  methods  is  inves cigated  by  comparing  results  with  oper¬ 
ating  boundaries  from  the  more  elaborate  analysis. 

In  this  volume  -che  differential  equations  of  motion  for  a  linearly  twisted 
rotor  blade  having  chordwise  mass  unbalance  and  operating  under  steady 
flight  conditions  are  derived.  The  motions  include  flapping  and  lagging  for 
the  articulated  blade,  as  well  as  flatwise,  edgewise, and  torsional  defor¬ 
mations  for  the  articulated  and  nonarticulated  blades.  The  fully  coupled 
aerodynamic  forcing  functions  are  based  on  quasi-steady  theory.  The  dif¬ 
ferential  equations  of  motion  are  expanded  in  terms  of  the  uncoupled  vibra¬ 
tory  modes  of  the  blade  in  order  to  facilitate  their  numerical  solution  on 
a  digital  computer. 
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FOREWORD 


The  work  presented  in  this  volume  is  part  of  an  effort  which  is  contained 
in  five  volumes.  The  work  was  performed  under  Contract  DA  1»U-177-AMC-332(T) 
with  the  U.  S.  Army  Aviation  Materiel  Laboratories,  Fort  Eustis,  Virginia. 
The  work  was  monitored  for  USAAVLAB3  by  Mr,  Joseph  McGarvey. 

This  volume  contains  a  presentation  of  work  performed  by  Mr.  Peter 
Arcidiacono,  of  the  United  Aircraft  Research  Laboratories.  The  resulting 
equations  of  motion  were  prepared  and  programmed  for  computer  solution  by 
Mr.  Russell  Berquist  of  Sikorsky  Aircraft.  The  method  of  computer  solution 
was  generated  independently  for  earlier  versions  of  the  equations  of  motion, 
and  was  extended  under  this  contract  to  include  non-coincident  blade  elastic 
axis  and  section  center  of  gravity  locations. 

Volume  II  presents  a  linearized  discreet  azimuth  classical  flutter  analysis 
for  rotor  blades,  with  an  appropriate  parameter  variation  study,  a  com¬ 
parison  with  test  data,  and  a  comparison  with  results  calculated  by  using 
the  method  of  Volume  I. 

Volume  III  describes  a  stall  flutter  analysis  based  on  the  calculation  of 
aerodynamic  work  during  a  cycle  of  blade  torsional  vibration,  using  two- 
dimensional  unsteady  airfoil  test  data.  The  analysis  was  used  to  generate 
stall  flutter  boundaries. 

Volume  IV  contains  the  results  of  a  study  of  static  torsional  divergence. 

A  set  of  design  charts  were  generated  and  the  effects  of  a  range  of  par¬ 
ameter  variations  are  presented.  The  results  of  the  static  divergence 
calculation  are  compared  with  results  calculated  by  using  the  method  of 
Volume  I. 

Volume  V  is  concerned  with  flapping,  flatwise  bending,  and  coupled  flap-lag 
instability.  A  single  degree  of  freedom  flapping  or  flatwise  bending 
analysis  was  used  to  investigate  a  wide  range  of  parameters.  Comparisons 
were  made  with  the  results  of  the  more  elaborate  method  of  Volume  I .  The 
method  of  Volume  I  was  also  used  to  determine  the  coupled  flap-lag  response 
of  a  rotor  to  a  number  of  sudden  control  changes. 
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INTRODUCTION 


The  development  of  analyses  for  the  prediction  of  rotor  performance,  blade 
motions,  stresses  and  loads  has  been  progressing  for  many  years.  The 
equations  of  motion  derived  in  this  part  represent  an  important  additional 
step  in  the  refinement  and  generalization  of  these  methods. 

The  analysis  used  to  develop  the  performance  charts  of  Reference  1  deter¬ 
mines  rigid  blade  flapping  motion  by  a  step-by-step  timewise  integration 
of  the  flapping  equation  of  motion.  The  solution  finally  converges  to  a 
cyclic  pattern  if  a  steady-state  condition  is  being  analyzed.  The  analysis 
is  fully  capable  of  handling  a  transient  condition,  however,  since  the 
equations  of  motion  are  solved  with  arbitrary  starting  values. 

Following  the  development  of  a  rigid  blade  analysis,  continued  improvements 
in  helicopter  performance  objectives  and  the  need  for  prediction  of  flex¬ 
ible  blade  loads  and  stresses  made  the  development  of  a  solution  for  flex¬ 
ible  blade  motions  necessary.  The  result  of  this  development  utilized 
rotating  blade  natural  vibration  modes  as  elastic  degrees  of  freedom.  The 
use  of  these  orthogonal  or  "normal"  modes  gave  rise  to  the  designation 
"Normal  Mode  Transient  Analysis".  The  analysis  is  similar  to  the  rigid 
blade  analysis  mentioned  above,  in  that  it  is  essentially  a  so-called  start¬ 
ing  value  problem,  in  which  the  differential  equations  of  motion  are  inte¬ 
grated  on  the  basis  of  some  set  of  starting  boundary  values.  When  a 
steady-state  condition  is_jbeing  analyzed  the  integration  proceeds  in  small 

but  finite  timewise  steps;  after  a  number  of  rotor  revolutions,  the  pre- _ 

dieted  motions  will  become  cyclic  within  a  desired  tolerance.  This  is  the 
usual  solution  desired,  and  performance,  load,  and  stress  calculations  are 
usually  based  on  these  cyclic  motions.  On  the  other  hand,  the  prediction 
of  rotor  behavior  following  an  arbitrary  initial  disturbance  is  a  basic 
capability. 

For  the  purposes  of  this  investigation,  the  above  analysis  was  extended  to 
provide  for  noncoincident  blade  elastic  and  center-of-gravity  axes.  The 
basic  differential  equations  are  supplied  in  detail,  in  order  to  define  and 
document  all  assumptions  completely.  The  provision  for  noncoincident  axes 
opens  entirely  new  areas  of  investigation  to  the  Normal  Mode  Transient 
Analysis . 


ASSUMPTIONS 


Nondimensional  quantities  are  used  extensively  throughout  this  report. 
Nondimensionalizing  factors  are  rotor  radius,  rotor  angular  velocity,  and 
mass  per  unit  span  at  some  blade  reference  station.  The  principal  assump¬ 
tions  which  were  made  in  the  analysis  presented  herein  are  listed  below. 

1.  The  aircraft  is  in  steady  flight  with  constant  rotor  angular 
velocity. 

2.  The  blade  has  an  elastic  axis  so  that  blade  deflections  can  be 
considered  as  the  superposition  of  two  orthogonal  translations 
of  and  a  rotation  about  the  elastic  axis. 

3.  Quasi-steady  aerodynamic  theory  is  applicable  with  the  exception 
that  apparent-mass  aerodynamic  effects  are  assumed  to  be  neglig¬ 
ible. 

h.  Radial  flow  effects  on  aerodynamic  forces  are  negligible. 

5.  Principal  blade  flexibility  effects  can  be  accounted  for  by  con¬ 
sidering  only  five  flatwise,  two  edgewise,  and  three  torsional 
vibratory  modes. 

6.  Blade  flap  and  lag  hinges  are  coincident  for  articulated  rotors. 

7.  Tlie  local  center  of  gravity  is  assumed  to  lie  on  the  major 
principal  axis  of  the  section. 

8.  The  blade  is  linearly  twisted  along  its  span. 

9*  The  following  quantities  can  be  assumed  to  be  small  in 
comparison  to  unity: 

a.  Flap  and  lead  angles  (in  radians)  and  their  derivatives. 

b.  Ratio  of  elastic  deflections  to  rotor  radius  and  their 
derivatives . 

c.  Ratios  of  chordwise  distances  (i.e. ,  chord,  center-of- 
gravity  offset,  etc.)  to  rotor  radius. 

d.  Built-in  twist  (in  radians). 

e.  Ratio  of  flap-lag  hinge  radial  distance  from  a  center 
of  rotation  to  rotor  radius. 

f.  Reciprocal  of  Froude  number  based  on  rotor  radius 

(  g/flzR  ). 

g.  Ratios  of  blade  thickness  dimensions  to  chord. 
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10.  On  the  basis  of  assumptions  8  and  9,  the  following  types  of  terms 
in  the  equations  noted  can  be  neglected  as  higher  order: 

a.  Flatwise  and  edgewise  bending  equations: 

(1)  Second-order  products  of  elastic  coordinates. 

(2)  Third-order  products  of  elastic  coordinates, 
chordwise  distances,  flap  angle,  lead  angle, 
built-in  twist,  and  flap-lag  hinge  offset. 

b.  Torsional  equation: 

(1)  Third-order  products  of  elastic  coordinates. 

(2)  Fourth-order  products  involving  elastic  coordinates, 
chordwise  distances,  flap  angle,  lead  angle,  built- 
in  twist,  and  flap-lag  hinge  offset. 

c .  Flap  angle  and  lead  angle  equations : 

Second-order  terms  involving  products  of  elastic 
coordinates,  chordwise  distances,  and  built-in 
twist. 

d.  Section  velocity  equations : 

(1)  Second-order  products  of  elastic  coordinates. 

(2)  Third-order  products  involving  the  elastic  coor¬ 
dinates,  chordwise  distances,  flap  angle,  lead 
angle,  built-in  tvist,  and  flap-lag  hinge  offset 
as  factors. 

e.  All  equations: 

The  spanwise  component  of  acceleration  due  to 
gravity. 
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AXIS  SYSTEMS 


The  axis  systems  employed  in  the  analysis  are  shown  schematically  in 
Figure  1  and  are  discussed  "below. 

NONRCTATING  ROTOR  HUB  AXIS  SYSTEM:  X,.V,.Z. 

This  axis  system  has  its  origin  at  the  rotor  huh  center  and  is  nonrotating 
with  .respect  to  the  aircraft.  Z,  lies  along  the  rotor  shaft,  positive 
up;  x,  is  normal  to  the  rotor  shaft,  lies  in  the  plane  formed  by  Zj  and 
the  remote  air  velocity  vector,  and  is  positive  aft;  y,  is  orthogonal  to 
X(  and  Zj  and  is  positive  toward  the  advancing  side  of  the  rotor.  This 
axis  system  is  Newtonian  in  nature, having  at  most  a  steady  translational 
velocity. 


ROTATING  ROTOR  HUB  AXIS  SYSTEM:  x2 ,y2.  ,z2 

This  axis  system  has  its  origin  at  the  rotor  hub  center  and  rotates  with 
the  rotor.  Its  coordinates  are  related  to  those  of  the  "1"  axis  system 
above  by  the  following  equations: 


xl  = 

x2cosiif/  —  y2  sini// 

(1) 

Vi  = 

y2  cos^  +  x2  sirn^ 

(2) 

2I  = 

z2 

(3) 

OFFSET  AXIS  SYSTEM:  x3 .V3  .z3 

This  axis  system  has  its  origin  at  the  coincident  flap 
radial  position.  Its  coordinates  are  related  to  those 
system  by  the  following  equations: 

and  lead-lag  hinge 
of  the  "2"  axis 

X2  = 

=  x3  H-  e 

(4) 

1 

v2 

=  V3 

(5) 

z9 

=  zx 

(6) 
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BLADE  LEAD  ANGLE  AXIS  SYSTEM:  *4  ,y4  ^ 

This  axis  system  rotates  about  the  Zj  axis  with  an  angular  velocity 
dS/dt  .  Its  coordinates  are  related  to  those  of  the  "3"  axis  system  by 
the  following  equations : 


X3  =  x4cos8  -  y4sinS 


(7) 


y3  =  y4  cos8  +  x4sinS 


(8) 


z3  =  z< 


(9) 


Note  that  8  is  equal  to  the  built-in  lead  angle  for  nonarticulated 
rotors , 


RIGID  BLADE  AXIS  SYSTEM:  X*  ,y5  ,Z5 

This  axis  systei  rotates  about  the  y4  axis  with  an  angular  velocity 
d/3/dt  .  Its  coordinates  are  related  to  those  of  the  "H"  axis  system  by 
the  following  equations: 


x4  =  x5  cos/3  -  z5sin/3 


(10) 


V4 


*5 


(11) 


z5  cos/3  +  x5sin/3 


(12) 


Note  that  the  flap  angle  <3  is  equal  to  the  built-in  coning  angle  for 
nonarticulated  rotors. 

RIGID  BLADE  SECTION  AXIS  SYSTEM:  x6  ,  y6  ,Z6 

This  section  axis  system  is  obtained  by  a  translation  along  *5  through 

r*  r  ow/1  o  nn  oVnnt  X-  +VirAiirrV»  1  HOfll  nffld  Til  t.oh 

U  UJ.O  VUUUV>  I  *  V  VV»VAW*«  MMVMW  5  w»»«  w  -  —  W - -  X- - 

angle  Q  ,  where  Q  includes  the  pitch  angle  resulting  from  control 
system  inputs  as  well  as  that  resulting  from  built-in  twist.  The  "6"  axis 
coordinates  are  related  to  those  of  the  "5"  axis  system  by  the  following 
equations : 


x5  =  x6  +  r 

(13) 

y5  =  ys  cos  9  -  z5sin 0 

(111) 

z5  =  z6cos  0  +  y6  sin 0 

(15) 

BENDING  TRANSLATION  AXIS  SYSTEM:  X  7  .  v  7  .z  7 

This  axis  system  is  obtained  by  bending  translations  of  the  elastic  axis 
in  the  y6  and  Zg  directions  (section  principal  axis  directions  before 
clastic  deformation) .  The  "7"  axis  coordinates  are  related  to  the  "6" 
axis  coordinates  by  the  following  equations: 

“6  5  *7 

(16) 

u 

+ 

< 

1% 

(17) 

z6  =  z7  +  we 

(18) 

Eq.  (16)  implies  that  the  change  in  length  of  the  blade  due  to  elastic 
displacements  is  negligible. 

-BENDING  ROTATION  AXIS  SYSTEM:  Xg  ,yg ,zg 

This  axis  system  is  obtained  from  two  small  rotations  which  occur  when  the 
blade  undergoes  the  elastic  be'  ding  translations ,  We  and  ve  .  The  rota¬ 
tions  are  assumed  to  be  small,  so  that  the  order  in  which  they  occur  is  of 
no  consequence.  The  "9"  axis  system  coordinates  are  related  to  those  of 
the  "7"  axis  system  by  the  following  equations: 

x?  =  Xg  -  z9  Xj  —  ygX2  (19) 


l-f  -  Zg  +  XgX,  (20) 
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(21) 


VELOCITY  AND  ACCELERATION  VECTORS  IN  THL  RIGID  BLADE  AXIS  SYSTEM 


To  develop  the  blade  equations  of  motion,  the  velocity  and  acceleration 
vectors  for  a  blade  mass  element  in  an  appropriate  reference  axis  system 
are  required.  A  convenient  reference  system  and  one  used  herein  is  the 
Rigid  Blade  Axis  System  ("5"  system).  Inasmuch  as  this  is  a  reference 
frame  moving  relative  to  Newtonian  space,  care  must  be  exercised  to  account 
for  all  velocity  and  acceleration  components.  The  approach  employed  herein 
assumes  that  the  linear  and  angular  velocity  vectors  of  the  Nonrotating 
Rotor  Hub  Axis  System  ("1"  system)  are  known  along  with  the  gravity  vector 
(from  specification  of  the  trim  flight  condition  of  the  aix'craft).  The 
velocity  and  acceleration  vectors  in  any  other  axis  system  can  then  be 
derived  by  classical  vector  techniques  as  discussed,  for  example,  in 
Reference  2  (Article  12.3).  The  velocity  and  acceleration  of  any  point 
defined  by  coordinates  measured  in  a  frame  of  reference  which  is  both 
translating  and  rotating  can  be  shown  to  be 


v  =  Vq  +  +  w  X  P  (27) 


g  is  the  gravity  acceleration  vector 

A  ^ 

i ,  j , k  are  the  usual  unit  vectors 
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Eqs.  (2T)  and  (28)  can  be  expanded,  if  the  above  definitions  are  used,  to 
yield  the  following  general  expressions  for  y  and  a  : 

V  =  +  x  +  uy  -  Wzy)  +  j{\^  +  y  +  wzx  -  w*Z) 

0 

+  k(vo  +  z  +  (jj, ty  —  wys) 


-a.  To  00  o  o 

=  ‘  LV0X  wzv0y+  wyvOj+  x  ~  2wzy  +  (wy  +  wxa/j)z 

+  2w^ +(-wz+ wxwy)y-(w/+w|)x  +  gx  ] 

r.  00  0  .0  . 

+  j  w,V(),+  Wz v0i^"  '  ”  2wrZ  +  K  +  wywx)x 

•L  y  z  22  "l 

+  2wzX+(-wx+  aJyW2}Z~{it>z+uv)y  +  Qy  j 

r  ijo  o  q 

+  k  [\Lj-  wyv0^  +  <jx  vQ^+  z  -  2  cuyx  -r  (tox  +  U)z  wy)y 
+  2  a>*y  -H-SyJ-  i«»Ia)x)x-{u>x2+  wf)z  +  gz  j 


where  (°)  =  -iil-i 
dt 

Eqs.  (29)  and  (30)  apply  to  any  axis  system.  For  application  to  the  Rigid 
Blade  Axis  System  all  quantities  should  be  subscripted  5>  and  expressions 
for  v  y  u>  y  x  >  v  >  an<^  Z  are  required.  These  are  derived 
below.  °5  5  5  *5  5 

The  velocity  vectors  o5  and  oj5  can  be  obtained  from  the  correspond¬ 
ing  vectors  for  the  "1"  axis  system  (  v  and  )  which,  for  steady 
flight,  are  1  1 


v  =  v  j  f  v  k 

vo,  o¥  1  o,  1 


a/,  =0 
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When  Eqs .  (31)  and  (32),  and  the  coordinate  relations  of  Eqs.  (l)  through 
(12)  are  employed,  and  when  Eq.  (29)  is  successively  applied,  the  follow¬ 
ing  expressions  for  vQ  and  £  result: 

5  5 

2  2 

%  ='fl{voil,cos*(|- T  ■  "f")  '  +  +  £v0J 

+  '5  {  _voX|(,_  f-)sil>  *  4  efl-v0<Scosi(i} 

+  k5  lv0Zi(*“  T ) -VoX|£cos*+VOx(/3  8sin*}  (33) 


cu5  =  i5  [/?(&+§)]  -4T5+  [a(.-f)  +s]T5 


J.f  Eqs.  (13)  through  (26)  are  used,  the  coordinates  X.  ,  y.  and  7 
can  he  expressed  in  tenns  of  the  blade  section  coordinates  (  X.n  y  5 

,°  )  as  f ol Lows  (assuming  blade  sections  having  infinitesimal  spanwise 
thickness  so  that  x  =0  )• 


r  ~  zio^i  ~  y. 


rIOX2 


(35) 


Vs  =  (ve+yio-zio0e)cos0-  (we+z,o+y|o0e)sin6 


(36) 


z5  -  ( we  T  zio  +  cos6  +  (ve  +  y|Q  -  zloee)  sin  6  (37) 


Substitution  of  Eqs.  (33)  through  (37)  into  Eq.  (30)  yields  the  following 
expression  for  the  acceleration  vector  in  the  "5"  axis  system: 

°5=  T5 {  -en2 -z,0  -  y10-^-  -2n[(ve-z,o8e)cos0-(ve+ylo 

-  Zo0e)  0 sin  8  -  (we+  ylo0e)  sin  9  -  (we+  z  10  +  yio0e)9cos 0  ] 

-  28|vecos  0-  (ve+y,o)03in0  -  wesin0  -(we+  z^cos©  ] 

+  (-js  +  /3ft2)  [(we+z1Q)cos  0+{ve+ylo)sin0]  -2/3[wecos0 

-  (wetz,o)0sin0  +  vesin  0  +  (ve+y,o)0cos0]  -s[(ve+ylo)cos0 

-  (we+z|0)sin  0]  -r[/32  +  n2  { I  -  /3 2 )  +  2&8+S2] 

-  il2[-Z|0-^r‘  ~  Yio  “dT  1  +  ^x5}  +  J5  {eil2S4  (Ove-Zio^)cos0 

-  2(ve  -z|o^e)0sin0-(ve+y|o-zio^)(0sin0+02cos0)  -  (we 

+  ylo°0e)sin0  - 2(we+yio0e)0cos0-(w5+z|o+y,o0e)(0cos0 -02sin0) 

-  2 (9  fi.  [ wfecos  0  -  (we  +  z^) 0  sin  0  +  vesin  0  +  (ve  +  yl0)  0 cos  0  ] 

,  0  ,0 

+  r(°S-  2/3/§il)  +  2.£l [- zl0^  - yl0-^r  ]  -•  2 pti  [( we+  z0) cos  0 
+  (ve+ylo)sin0]  -^[(Vg+y^-z^Jcos  0-(we+zlo+ylo0e)sin0] 

-  2il§[(ve  +  y|O)cos0  -  (we  +  z|o)sin0j+ gyjj  +  k5{/3&2e 
+  (we  +  y,o0e)cos0  -2{we+y|O0e)0sin0  -(we+zt0+  ylo0e)(0sin0+02cos0) 

0  +  (Vg-zio^)sin0  +  2(ve-z|o0e)0cos0  +  (ve  +  yio-z|o0e)(0cos0-02sin0) 

+  2/3il[vecos0-(ve+yjO)0sin0-wesin0  - (we+ zl0) 0cos  0 ] 

+  r  [/?•+  (B{Sl  +  28il) j  +  9zs  }  (38) 
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EQUATIONS  OF  MOTION  FOR  BLADE  ELASTIC  DEGREES  OF  FREEDOM 


Details  of  the  derivation  of  the  equations  of  motion  for  the  blade  elastic 
degrees  of  freedom  are  given  below.  Briefly,  the  equations  were  derived 
from  consideration  of  the  equilibrium  of  aerodynamic ,  dynamic ,  and  elastic 
moments  at  a  given  section  of  the  blade.  The  resulting  moment  equilibrium 
equations  were  differentiated  as  required  to  express  them  in  terms  of  the 
local  blade  loadings ;  the  elastic  deflections  cf  the  blade  were  replaced 
by  summations  of  appropriate  uncoupled  blade  vibratory  modes,  of  which  a 
finite  number  were  retained  in  the  final  analysis.  Stiffness  terms,  both 
structural  and  centrifugal ,  appearing  in  the  equations  were  replaced  by 
equivalent  natural  frequency  terms ,  and  the  orthogonality  properties  of  the 
assumed  modes  were  employed  where  possible  to  eliminate  modal  dynamic  cou¬ 
pling  terms. 

MOMENT  EQUILIBRIUM  EQUATIONS 

The  basic  moment  equilibrium  eouations  for  any  blade  section  are  of  the 
following  general  form: 


(Me)x|0  =  (Ma)x10  +  (39) 


(Me)yo  =  (MA)yD  +  (M0)y(o  (UO) 


(M«>zl0  =  (Ma)Zk)  +  (M0)Z)0  (141) 


where  Me  ,  MA  and  M0  represent  moments  due  to  elastic  defor¬ 
mation,  applied  aerodynamic  forces,  and  dynamic  body  forces  respectively. 
It  is  convenient  to  express  the  latter  two  moments  with  reference  to  the 
Rigid  Blade,  or  "5"  Axis  System.  When  the  axis  transformation  relation¬ 
ships  given  in  Eqs.  (l)  through  (2 6)  are  employed,  Eqs.  (39)  thro.igh  (hi) 
are  as  follows : 


0*2) 
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(43) 


(Me)y(o  =  My^  (cos  9-  9e  sin  9)  +  Mz^(sin  9  +  9e  cos  9 1  -X2  Mx 


(Me)Z|0=  MZ5(cos0-ee  Sin 0)  -  My5  (Sin  9  +  0e  cosS)-X(Mx  (UU) 


where  ^x5  ,  ^y5  ,and  ^zg  represent  the  sum  of  the  aerodynamic  and  dy¬ 
namic  moment  components  in  the  xg  ,  yg  ,  and  Zg  direction,  respec¬ 
tively  (e.g.,^xg  =  )x  +  ^0^X5  )•  Solving  for  from  Eq. 

(42),  substituting  the  result  into  Eqs.  (43)  and  (44),  neglecting  higher 
order  terms,  and  defining  a  total  pitch  angle 


©  =  0  +  9e  (45) 


yields  the  following  forms  of  Eqs.  (43)  and  (44): 

(Me)yio  =  My5  cos©  +  MjjSin©  -  ^Me)X|0  (46) 


cos®  -  Myjsin©  -  X,(Me) 


It  is  now  convenient  to  nondimensionalize  Eqs.  (42),  (46),and  (47)  by  using 
as  nondimens ionali zing  factors  the  rotor  radius,  R  ,  the  rotor  angular 
velocity,  £1  ,  and  a  reference  mass  per  unit  length,  rnQ  .  All  such 
nondimensionalized  quantities  will  hereafter  be  indicated  by  means  of  a  bar 
(e.g.,  7=  r/R,  Siy5  =  My  ■'  mQi22R3  ,  etc.);  further,  ^differentiation 

with  respect  to  nondimens ionll  time  will  be  indicated  by  (  x  ) ,  while  dif¬ 
ferentiation  with  respect  to  nondimens ional  radial  distance  will  be  indi¬ 
cated  by  (  )'.  If  these  conventions  are  employed  and  Eqs.  (22)  and  (23) 

are  used,  the  moment  equilibrium  equations  are  then 

(M  )  =  My  +(vg/-0,wJ(M  cos0  +  M  sin 0)+(wp  +  sV)(M7  cos  0 -&L  sin  0) 

e  x(0  xs  c  e  ys  ls  e  e  £s  ys 

(48) 
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My  COS  Q 
>5 


+  M,  sin  ©  -  (M  Jv 

Z9  e  XI0 


MJz  -- 

no 


Mz_cos©-My  sin© -(Me)x  (  we'  + 
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vte') 


(50) 


The  moments  associated  with  the  elastic  deformation  of  a  section  of  a  beam 
having  a  finite  rate  of  twist  along  its  span  and  a  center  of  tension  dis¬ 
placed  in  the  chordwise  direction  from  the  elastic  axis  have  been  derived 
in  Reference  3.  Expressing  the  results  of  Reference  3  in  terms  of  the  yve 
and  ve  elastic  coordinates  employed  in  this  analysis  and  neglecting 
higher  order  terms  in  accordance  with  assumptions  8  and  10  lead  to  the 
following  expressions  for  the  elastic  deformation  moments: 


(rae*yl0  =  “  ^  ^we  +  ^  ) 


(51) 


(Me)z10  =  EIz  (ve"-2weV) 


(52) 


(Me}x 
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GJ  0e  +  (B1  +  Be') 


df 


(53) 


The  last  term  in  Eq.  (52)  is  not  contained  in  the  results  of  Reference  3 
hut  rather  arises  when,  for  blade  balance  purposes,  a  counterweight  is 
employed  which  is  so  mounted  as  to  be  restrained  radially  by  the  blade’s 
tip  structure  and  laterally  by  the  blade's  local  structure.  This  type  of 
counterweight  is  then  constrained  to  bend  with  the  blade  and  is  in  com¬ 
pression  under  the  action  of  centrifugal  force. 

Combining  Eqs.  (51)  through  (53)  with  Eqs.  (1»8)  through  (50)  and  neglecting 
higher  order  terms  in  accordance  with  assumption  10,  results  in  the  final 
form  of  the  moment  equilibrium  equations  below: 


GJ0e'  +  (0,+ee)kA2/-  VC  =  +  ( ve' -  e'we)(My5  cose  +K/i2^ sinS) 

+(wa/+e,va)(M7  cose-fl„  sine) 

(5fc) 


V 


-Ely  (we"  +  2%‘e1)  =  My^cos©  +  M2sSin@-(Me)X|o(vt'  -  we$') 

(55) 

Ei2  (5j'-2we'8')-ea  jf\df  +  /V  5Xj  d{ 

°CW 

=  ii2jC0S@  -rayjsin®  -(Me)X|0(w'+ ve8') 

(56) 


The  last  two  terms  in  Eq3.  (55)  and  (56)  are  always  comparatively  small, 
and  are  neglected  hereafter. 

LOADING  EQUILIBRIUM  EQUATIONS 

t 

Flatwise  Loading  Equilibrium  Equation 

The  desired  flatwise  loading  equilibrium  equation  can  be  obtained  by  dif¬ 
ferentiating  Eq.  (55)  twice  with  respect  to  the  spanwise  coordinate  r  (if 
techniques  described  in  Reference  U,  pp.  352-353  are  used).  The  resulting 
equation  is 
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~[ETy  (we//+  2ve/0,)j//  =  ©"[-My^in®  +  MZscos®] 

+  (®')2  [~^y5C0S®  -  MZssin®j  +  2®^- My5sin® +Mz5cos®j 
+  (videos®  +  M2"sin®  (57) 


To  proceed  further,  expressions  for  the  "5"  axis  moment  components  are  re¬ 
quired.  These  are  derived  in  Appendix  I  along  vitn  their  derivatives  with 
respect  to  r  .  The  derivatives  of  ^y5  and  in  Eq.  (57)  are 

replaced  by  the  expressions  given  in  Eqs.  (110)  through °( ll4 )  of  Appendix 
I.  Then  by  employing  the  relations 


■5eaS,n® 


y5eaC0S©  A  - 


8'  +  v.' 


~  z5eaC0S®  =  -*•  +  (8')Z we  -2v‘8'-ve8" 

(60) 

y5e0COS® +f5easin®  =  veW-  ve<0')2  -  2wfiV-  we0" 

(61) 
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and  by  neglecting  higher  order  terras  in  accordance  with  assumption  10, 
Eq..  ( 57 )  becomes 


-  [ny-(Wew  +  2ve'0')f=  -  2©'{sin©jrrTSZ5d^+cos®^7TSy5d$}  -  vi£ 

-  SZscos©  +  5y5sin©  +  ^5(we'  +  \e)  +  2®cos0  [^yj' 

+  ©'yiOc9^x5,y10cg:o  +  ©"[-My5sin©+M25cos©] 

+  (©') 2[-My5C0S  ©  -  Mzs\n  ©]+  20/fsin  0(mA)y5-(mA)Z5  cos©] 

-  (m.)/  cos  ©  -  (m.)/  sin  © 

*5  *5 

(62) 


The  flatwise  loading  eauation  can  now  be  expressed  in  the  desired  form  if 
the  ©"  and  (@')2  terms  are  expanded  through  use  of  Eas.  (^5),  (55) 
and  (56),  and  if  higher  order  terms  are  neglected  in  accordance  with 
assumption  10. 


//  rT  jfy 

-[E!y(We/+  2veV)]  -  -2©' {sin  ©  JL  SZ5dC  +  cos©  J  Sy,.d£}  -  w/jL  S^df 


-  S^cos©  +  Sy&sin  ©  +  SXj.(we4vee/)  +20cos0[ffiky)^] 

+ ©  *oc< S0 )x5, y10  =0_6e eA^  Sx&dC  +  A \lTrdt 


+ 


2®/  [(mA^sin®  -  (mA)Zjcos©j  -(nriA^cos©  -(m^sm© 


(63) 
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Edgewise  Loading  Equilibrium  Equation 


The  desired  edgewise  loading  equilibrium  equation  is  obtained  by  differ¬ 
entiating  Eq.  (56)  twice  with  respect  to  r  .  The  resulting  equation  is 


[eiz  (ve  -  2weV)  ~eAf  SXgdC  +  AeA{Jl  df  ] 

°cw 


// 


=  [  M^sin®  +  My^cos©]  +  (Q^2  [-M^cos®  +  Myosin  ®] 

-  20/['M2(./sin©  +"MyJcos ©  ]  +  M2"cos  ©  -"Myosin® 


(6b) 


Substituting  (from  Appendix  I)  for  the  derivatives  of  My5  and  Mz  , 
using  Eqs .  (58)  to  (6l)  as  required,  and  neglecting  higher  order  terms5 
where  necessary,  allows  Eq.  (6b)  to  be  expressed  as 


[eIz(  V  2we0) -ej  \d(  +  A  «  ] 

r°cw 

=  -  ©"  [ M^sin  ©  +  My5cos  ©  ]  +  (©')2  [  -¥z.cos  ©  +  My5sin  ©  ] 

-  2©/{  -sin  ©  j_  TSy5d^  +  cos  ©  jf  \d?  }  +  \jf  \  d£  +  Sy5cos  © 


+  S7  sin  ©  - 


s  -  (-wee'+v^5+  [iioJSo^oj'-aesine^m)' 


+  2©  [  (mA)^sin  ©  +  ( mA)^  cos  ©  ]  -  (mA)z  cos  ©  +  (m  A)^  sin© 


(65) 
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Or,  expanding  the  ®;’  and  (®()2  terms  through  the  use  of  Eqs.  ( 1+5 ) , 
(55) »  and  (56)  and  neglecting  higher  order  terms  permits  Eq.  (65)  to  be 
placed- in  the  desired  form,  as  follows: 


=  -2©'  {  -  sin  ©  f.  \5d{  +  cos©  £  52£.dc}  +  ve^rTSX5dC+^y5cos© 
+  SZ5sin  ®-(-wee  + ve')SX5+  [y^So^y^J-^sin  B (kz*  m/ 


2®  [  (mA)  sin  ©  +  (mA)  cos  ©  j  -  (mA/  cos  ©  +  (mA) '  sin© 

*•  Za  y.  Cm  V* 


y5 
(66) 


Torsional  Loading  Equilibrium  Equation 

The  desired  torsional  loading  equation  is  obtained  by  differentiating 
Eq.  (5*0  once  with  respect  to  r  .  Prior  to  differentiation,  however,  it 
is  convenient  (l)  to  substitute  for  lWx  with  the  use  of  Eq.  (106)  of 
Appendix  I,  and  (2)  to  simplify  the  ^y5  5  and  Pz  terms  in  Eq.  ( 54) 
with  the  use  of  the  definitions  below:  5 


=  VeCosS-O'wgCos  0-w/sin  9-d\c  sin  8 


(67) 


75ea  =  vesin  8-  0wesin  8  +  wecos  8  +  8%  cos  8 


(68) 


The  resulting  expression  for  the  torsional  moment  equation  is  then 

GO  8.'+  («'+«.')  --  <m4y  sl5(y5m-  ygn) 

-  5y5(V?5»(F))}  d?  +  ^  V,oJ  <sD)Iscos  ©  -(s0)  sin  ©  j  d£ 

j  _  V 

-  jj  m(kZj2o+ky2|o)(©  +  ^)dC  +  ^T^-kZ|o)[^(i  +  2S)siri20 

+  {8e-0 )  cos20]d£  -  zj.  m  k^gsin  9  d£  -  2  J.  rnky^WeCos  8  d£ 


+  V*.'.+ V** 


(69) 
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Note  that 
functions 
indicated 
spect  to 
Eqs.  (55) 
torsional 


all  quantities  appearing  in  the  various  jntegrands  above  are 
of  the  dummy  variable  of  integration  £  ,  unless  explicitly 

to  be  a  function  of  r  .  Differentiating  Eq.#/(69,  with  re- 
r  and  using  the  definitions  of  y5"  and  2.  along  with 
and  (56)  to  eliminate  Mz  and  My  80  results  in  the  final 
loading  equation  below:  5  5 


Jyiocg:° 


[Gj8e'+(8'«;)k lj_  ’sX5dj -(mA)Xs  -y,0J(S0)2  “S®-  (SD)y!>sm  ®] 

?  p  XX  X  — -  ?  —  2  f  I  x  X  1 

+  m(kz  +  kyl(©+/9)-m(ky)o-kZ|o)[^  (l  +28) sin  28  +  (8e-/3)cos28  J 

+  2mk^Qsin8(Ve/-  8we')  +  2mky%os8(We  +  9%)  + 

+  (we//+2ve/8l [-eA  f  +  AeAcwi^  S^df] 

*  °CW 

+  ( EIZ  -  EIy ) ( we"  V"  2  '  +  2ve"ve/8/) 


(TO) 


MODAL  EQUATIONS  OF  MOTION 

The  nreceding  sections  of  this  report  have  been  directed  toward  the  der¬ 
ivation  of  the  fundamental  differential  equations  of  motion  (Eqs.  (63), 

(66),  anc  (70)) governing  the  forced  response  of  the  elastic  motion  of  a 
twisted  rotating  beam.  These  equations  have  generally  been  derived  in 
terms  of  the  local  dynamic  and  aerodynamic  shear  forces  and  moments.  Ex¬ 
pressions  for  the  dynamic  shear  forces  are  given  in  Eqs.  (102),  (103)  and 
(108)  of  Appendix  I,  while  the  aerodynamic  shear  forces  and .mechanical 
damping  forces  are  derived  in  Appendixes  II  and  III  respectively.  Examin¬ 
ation  of  the  dynamic  and  aerodynamic  shear  force  expressions  indicate  that 
closed-form  solutions  to  the  blade  equations  of  motion  are  not  possible 
because  of  the  nonlinearities  present.  Therefore  solutions  must  be  ob¬ 
tained  by  approximate  means.  The  particular  approach  employed  herein  is  a 
modal  approach  wherein  the  elastic  coordinates  we  ,  ve  ,  and  ^  .  , 
are  expressed  as  finite  series  summations  of  assumed  radial  shape  functions 
(mode  shapes)  with  each  suitably  scaled  by  time-dependent  generalized  coor¬ 
dinates.  The  series  expansions  employed  are 


=  I  yWn(r  )qWn  («^) 


n=l 


(71) 


ve(Tv/>)  =  iyVm(r)qVm  (*} 

m=| 


(72) 
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(73) 


3 

0e<7W')=  I 


Note  that  the  number  of  terms  retained  in  the  expansions  of  Wg  ,  vfi  , 
and  0e  are  5,  2,  and  3,  respectively.  The  radial  shape  functions 
(hereinafter  referred  to  as  mode  shapes)  employed  in  this  analysis  cor¬ 
respond  to  the  natural  vibratory  mode  shapes  of  an  untwisted  blade  mounted 
on  an  unaccelerated  rotor  hub  and  operating  at  zero  pitch,  flap,  and  lead 
angles  and  at  a  rotor  angular  velocity,  ft  .  As  indicated  in  the  anal¬ 
ysis  which  follows,  the  substitutions  of  Eqs.  (71)  through  (73)  into  Eqs. 
(63),  (66),  and  (70)  permit  these  equations  of  motion  to  be  expanded  into 
a  system  of  simultaneous  differential  equations  with  the  generalized  co¬ 
ordinates  as  the  unknowns.  Further,  because  of  the  orthogonality  prop¬ 
erties  of  the  particular  mode  shapes  employed,  the  resulting  system  of 
equations  can,  to  a  large  extent,  be  dynamically  uncoupled  as  far  as  the 
blade  elastic  degrees  of  freedom  are  concerned.  As  a  result,  solution  of 
the  equations  through  numerical  integration  techniques  is  greatly  facil¬ 
itated.  Details  of  the  procedures  followed  are  given  below. 

Flatwise  Modal  Equation 

Substituting  Eqs._(7l)  through  (73)  into  Eq.  (63),  setting  ^mAV5  , 
and  (Sft)x  equal  to  zero  in  accordance  with  the  results  of 
Appendix  II,  and  employing  the  uncoupled  free  vibration  form  of  Eq.  (63), 
i.e. , 


>  El 
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•)//  _ 

+  m 


(e  +  r) 


(7*») 


2? 


permits  Eq.  (63)  to  be  expressed  in  the  following  form: 

yw„  I*"*.*  -2 [si, q,„  ]  =  -2©'{sin®  l Ts  df 

111  f  5 

+  cos®  f_  Sy5  [  Jjs0  l<5df  -f_  m(e+f)df] 

-  ,,  e'l 

5  — 

x  _  T. T  _ 

26cos  9  (m  ky(Q  }  +  0y'OcglsO  )x5,yl0  =o_0e  eA  J  (So)x5d6 


«w  */ U  Q  y 

5  rn  vfn  vm 


r  •— 


+  A%.  <L  'Soi„dT 


cw  ~  *7  m 
°cw 


(75) 


Note  that  in  Eq.  (75)  and  in  succeeding  equations,  a  shorthand  summation 
notation  is  used  whereby  the  summation  limits  are  not  explicitly  indicated 
but  rather  are  understood  to  be  those  of  Eqs.  (71)  through  (73). 

The  integrals  containing  the  Z5  and  yg  shears  can  be  simplified  con¬ 
siderably  by  noting  that  when  the  response  of  a  beam  is  represented  by  a 
summation  of  modes ,  the  shears  at  any  point  on  the  beam  may  be  expressed 
in  terms  of  the  modal  inertial  loadings  and,  hence,  in  terms  of  the  modal 
amplitudes  and  associated  modal  natural  frequencies  (Reference  5,  p.  64l) . 
Thus ,  the  following  relations  are  valid  for  a  limited  mode  analysis : 

/!_  _  __  «T _  _  Jj_  _ 

l  %<*£=  ?wwfqwnXmywn  cos0d£+£  (1+  cuv5qVm/myVmsin9d£ 

r  r  m  r 

(76) 


5in0d£  +  jT  (l+cuv2)qVm  /TmXVmcos0de 


(77) 
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The  I’s  in  the  second  summations  above  result  from  the  centrifugal  shears 
associated  with  the  uncoupled  edgewise  modal  displacement.  Employing  Eqs. 
(76)  and  (77)  and  the  definitions  of  5,  ,  §„  ,  (§o)x5  , 

and  ^  from  Appendixes  I  and5II  permits  Eq.  (75)  to  be 
expressed  in  the  following  form  (after  higher  order  terms  have  been  elim¬ 
inated)  : 


"  y"nq*n  Wwn  +  2  fol*.  QV(n^,  ]  =  20,{  £  0)*^  J_  m 
+  ?yelllr)qek+eB(r,)C^+  ^(l+^v2  ) Qv  /Tmyv  d^} 

K  *  *  m  m  m*i  vm  i 

+  (-9b+«b(7>!  « +£< 1  } 

+  5ywX,  ^Im[2?5-2yIOclosin8o]d{-SIrw'nqwj2‘rS-27,Otje0sin80] 

-  m 8, 1  r,„ qVm-280coseo(  Sk>|*),-(9|+  Irs'  q^ ) 

+  ( VVW:5  cos  ®  -  ( sA  -  $L[)ygsin  ®  ~  ^  £  cos  0Q- e  8  sin 0O 

p  Z  * 

+?ywn[qXwn+qwn(-4  "2/30+0  sin20o-i<l  -  cos 20+  20Bsir.20o> 

-8  <1-  cos20o»]+ly|Oc +ZyVtnqVm)[/3(  i-COS20o)-^(sin20o  +  20ecos20o) 

+  Ssin20o  +  0J+  ZyVfn QV(n  [ii0+2fi]  +  IyVm\2y.OcSin0o 
+  £  y0k  ^biOcg005  26o  ~  r  ( /3+/3)  sin0o-  rScos^  +  y^Iy^q^ 

yel V  [‘6a|  +  A\w  /_  d  c]  +  r  ( /3+/3 )  ( cos  0O-  0^in 0o> 


*  "  — —  __  _  XX  ^  ^ 

+  2/38  r  cos90+ qzcqs90-  gy5Sin0o-r8(sin0o+8Bcos0J  +  2r/3/9  sin0o} 
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The  equations  of  motion  for  each  flatwise  mode  can  be  obtained  by  consid¬ 
ering  the  work  done  during  a  virtual  displacement  of  each  mode,  rormaiiy, 
this  invloves  multiplication  of  Eq.  (78)  by  >Wj  ?  integration  over  the 
blade  span,  and  application  of  the  orthogonality  property  of  the  assumed 
uncoupled  natural  modes ,  namely , 


(79) 


This  procedure  yields  the  following  equation  of  motion  for  the  i^*1  flatwise 
mode : 


0  = 


/  >w,  [(SA)Z5c°s©-(SA)y5sin©]  d r - MpR cos 0P R(yw' )_sq  +  si  n0LOMLD(yw'L;o 

-  cos0FDMFO(yw^io-C|.  [qWj  +  qWl  (ww2-0o2-2/30o  +  /3  sin 20o-  -i-  (  i 

+  28)  (i  -  cos20o))]  +Zqw  [-20^^09  +0  sin20<£  +  28(C  -C  ) 

v  v  v  i y  1  wy/  7  4 1  j/ 

'  20osir8o(C66-c69  )i,i +  £  qvc(c3i)[-S0-^  (l-cos280;--|-(i+2S)sin291j 

-  8, COS  280(C5-O.T5C3)i/  C27i  o+  28,  ( ,  +^*)  ce  ^-8,0,  J 

5J-2SineoC64,,-2(«o^!C3J  +*rrq,| 


"?%,C62.  ,  +  ?%,[C: 


j  “i  62 1, i  j 


57'  -C62,||COS20o  +  C56,;j 


(( /3  +/3  )  sin0o+8cos0o) 

~  C®*ii j]  ~CiO|[e  (/3cos0o~8sin0o)  +  ^cosOo-g^sm©,,]  -q2j[(/3 
+  ^+2/38)  cos 0O  +  ( 2/3/3  -8 )  s  in0o]  +  C|3 ( [  ( 0  +/3 )  si  n0Q  +8cos0o] 

-C63j[0o+  /3  { I  -COS20O)  +  ( I  +  28  )  Sin  20o]  -0,  COS20d(C6oi-O.75C63|)-0|C6oj 

(80) 


where  the  various  modal  integration  constants  (e.g.,  C(.  )  are  defined 

in  Appendix  IV.  Eq.  (80)  above  is  the  final  flatwise  modal  equation  of 
motion.  Note  that  in  this  equation,  as  in  the  final  edgewise  and  torsional 
modal  equations  to  follow  in  the  next  sections ,  the  modal  subscript  nota¬ 
tion  has,  for  convenience,  been  standardized  so  that  the  i  subscript 
refers  to  flatwise  modal  quantities,  the  p  subscript  to  edgewise  modal 
quantities,  and  the  j  subscript  to  torsional  modal  quantities.  Where 
products  of  flatwise  modes  appear,  j'  subscripts  are  introduced  to  dif¬ 
ferentiate  between  terms  associated  with  the  various  summations.  The  p' 
and  j'  subscripts  serve  similar  purposes  for  the  edgewise  and  torsional 
modes.  Furthermore,  Eqs.  (155) >  (156) ,  and  (l6l)  have  been  substituted 
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V 


for  the  flap  damper,  lag  damper,  and  pushrod  forces,  respectively;  and  the 
integrations  over  7  have  been  performed  to  obtain  the  corresponding 
moment  contributions  defined  below: 


“lo"  ( s+  £  (x''»,,So,v1>cosSl0" q*P8osinS'-i>) _ 

?  ^W|  S'n  ®LD*" ^W|®cP°®®L0^  ] 


(81) 


^FD  "  m0Il0R3  [  ^  +  £  ^wiV.o  ^wi  C0S  ^FO  “  ^Wi  ^o  S’n  ^FD  ^ 


1  (y') 

p  VP  F  = 


X 

+  e0  cos 
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Edgewise  Modal  Equation 

The  desired  edgewise  modal  equation  is  derived  in  a  manner  entirely  anal¬ 
ogous  to  that  followed  in  deriving  the  flatwise  modal  equation.  The  free 
vibration  identity  and  the  orthogonality  condition  for  the  uncoupled  edge¬ 
wise  modes  used  in  the  derivation  are 
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myvm^vm2 
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J T  i»  rT  _ 

1 1  n>rVmyVpdr  -  l  mx^dr 


(85) 


The  resulting  edgewise  modal  equation  is  (with  the  standardized  j  ,  j 
p  modal  subscripts) 

o  =  jf  yv  [(SA)y3  cos©  +  (§A)Z5Sin©J  dr  -(yVp)Ft0  MPR  sin0PR 

-(/')  sin0FDMFO  -  (yVpl  cos0FD  Ml0  -  Cl4p  [qVp  +  qVp  (wv2 

VP  r=0  r  r=0 
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—  a  2 


#0  ~  2(380  -  0sin20o  -  y(l  +  2  8)(i*cos20o)  +■  1^  j 
+  j?  qv["C34PiP/  sin2  0o  +28  (C20-C|8)p  ^-2^0sma0(C77+C5l+C^Pip/ 
+  20,  C2(p  p/  (|  +  o5Vp2)]  +  2,  V  2  COS  0Q  (C5(+C96t2C79)p)P/ 

-2^1  q^qwt^vf  C75liP,j  +  J  Qw,  [c»s||P-  2e,aJwf  C22)iP 

—C3 i  p{d0-k  (1+  cos2 90)+y  (l  +  28)sin2  0o^  +  0,C5. 
0,COS20o(C5-O.75C3).p-20oCOS0o(C53+C9S+C8|  )(p  ] 


5i,P 


-  -  ,|P  j 

T<-  ^wl  C3.  +  ^  "  S'n^C53+C99+C8l  l  ]+2#0XqflC74  . 

I  1  'ip  UP  j  i  P,J 

+  ? q0,[‘C74Pifin2®o“cioop/  ^+^)cos0o  -Ssin0o)]-ci5Je(/3sin^ 
+  8  COS0Q)  +  gZs  sin0o+  gyg  cos0o]  -  +  2  08)  sin 0O 

+  (8-2/30)  COS0J  -C33p[(0+0)cos0o-  8sin0j+  c7oj0o+200o 
+  *  sin26o  +  T  (l  +2S)d+cos20o  )]-sin20oC76p+e(c5O+c72 
+  C23^p  +(  l+28  )  (C 52+  C73+C TO+  C24) p -  20o  sin0oC9lp  (eg) 


The  modal  integration  constants  are  defined  in  Appendix  IV. 

Torsional  Modal  Equation 

The  torsional  modal  equation  is  derived  as  follows:  First ,  Eqs.  (71) 
through  (73),  (l**7),  and  (1^8)  are  substituted  into  Eq.  (70)  and  ky2Q 
is  neglected  compared  to  kz(2  (a  valid  assumption  for  blades  having 
thin  ai^ foil- type  cross  sections).  Next,  the  torsional  free  vibration 
identity  for  uncoupled  vibration, 


V 


and  the  shear  force  definitions  of  Appendix  I  are  substituted  into  the 
resulting  equation.  The  final  torsional  modal  equation  is  then  obtained  by 
multiplying  the  torsional  loading  equation  (as  modified  by  the  above  oper¬ 
ations)  by  &Q.  ,  by  integrating  over  the  blade  in  the  7  direction, 

and  by  applying Hhe  orthogonality  condition  for  the  assumed  uncoupled 
torsional  vibratory  modes,  namely, 
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The  resulting  torsional  modal  equation  is  (with  standardized  modal  sub¬ 
script  notation  being  used) 


0  = 


*  KL  dr+C,,+C0S28C8o-2q  C  -C  qe 

i  i  i  i,l  1  i 


t 

+  C36jq0j  (|-a'0j2  -cos20o)-C44  [fl0+  T  { 1+28  )sin20o+£  (l-cos20o)] 

“  S  qwi'C38i,i(i+0'^,qvpqvp'C39plP<j 

+  f  qw,  [2fl°SI  n  0°  C 110 1 ,  i "  C86jt  f  2  8  ^ j"  ® C83 i  (  -  (I  +  2 8 )  C84|  - c62 . ^  ,<-0o  2 

"2^0O+  P  sin20Q  -  —■  (l  +  2S)(l-cos20o)>  +  sin20oC95  ] 

',iJ 

/2C74„ J  4+0  )]+Iq„J-29|  C87nT8,  <^rC„.  ,<« 

X 


P.I 

X 


"P.J  '^p,l 


XX 

0 


.  112,  ,/ ^Wj^/S'n^o 


+/3  ( i  -cos20o)  +  -j-  ( i  +  2  8 )  sin  2 6^  -  cos  20OC„4  | 

p,iJ 

-ZZ  2C„,  qw  (q  cos0o  -qv  0osin0o)+H2C1 

i  p  *»p,i  ■  p  p  i  j/ 

-qw,A  C°se0)-c92j  [i(0cos0o-  8sin0o)  +  gZ5cOs0o-ty5sin0Q| 
“c93,  [(^+/3+2/38)cos0o-(  8-2/3/3)sin0o]  +Q*  [  (/3+£)sin0, 
+Scos  0O]+Iq0  C„  [(|+^)sinfl0- 8cos£j 

j/  J  Iff  J 


(89) 


l 


28 


EQUATIONS  OF  MOTION  FOR  BLADE  FLAP  ANGLE 
AND  IEAD  ANGLE  DEGREES  OF  FREEDOM 


FLAP  ANGLE  EQUATION 

The  flap  angle  equation  can  be  obtained  directly  from  Eq.  (98)  by  setting 
the  r  integration  limit  to  zero  to  obtain  the  moment  about  the  flap 
hinge.  This,  in  turn,  is  equated  to  zero.  With  the  elimination  of  r 
items  in  the  integrand,  £  (the  dummy  variable  of  integration  in  the  r 
direction)  can  be  replaced  by  r  to  yield  the  following  equation: 
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0  1  /  {<"Wy5  '  V  + +  Vioc9  sin®  (s„  ) 


x5'yioc,=° 
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Substituting  the  shear  force  and  moment  expressions  from  Appendixes  I 
through  III,  and  the  modal  expansion  equations  (Eqs.  (71)  through  (73)) » 
neglecting  higher  order  terms  in  accordance  with  previously  defined  assump¬ 
tions,  and  noting  that  the  uncoupled  flatwise  and  edgewise  mode  shapes  are, 
for  small  hinge  offset,  essentially  orthogonal  with  the  function  r  , 
results  in  the  following  desired  flap  angle  equation  of  motion: 
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SUMMARY  OF  PRINCIPAL  EQUATIONS 


In  view  of  the  length  of  the  analysis  described  in  this  report,  a  summary 
of  the  principal  equations  is  given  here. 

The  blade  equations  of  motion  are  as  follows : 

Flap  Angle  -  Eq.  (91) 

Lead  Angle  -  Eq.  (96) 

i^h  Flatwise  Bending  Mode  -  Eq.  (30) 

p^*1  Edgewise  Bending  Mode  -  Eq.  (86) 

jth  Torsional  Mode  -  Eq.  (89) 

The  above  equations  of  motion  are  given  in  terms  of  section  aerodynamic 
forces  and  moments,  flap  damper,  lag  damper  and  pushrod  moments,  and  modal 
natural  frequencies  and  integration  constants.  The  section  aerodynamic 
forces  and  moments  are  given  by  Eqs.  (l^l)..  (1^2),  and  (lU6);  the  damper 
and  pushrod  moments  are  given  by  Eqs.  (8lJ  through  (83);  and  the  modal 
integration  constants  are  defined  in  Appendix  IV.  The  required  uncoup] 2d 
blade  mode  shapes  and  frequencies  can  be  determined  from  an  appropriate 
eigenvalue-eigenvector  analysis. 
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CONCLUDING  REMARKS 


1 

The  differential  equations  outlined  in  this  volume  are  an  extension  of  an 
existing  Sikorsky  Aircraft  advanced  method  for  the  prediction  of  rotor 
loads,  stresses,  and  performance.  In  most  cases,  these  predictions  will 
he  good  approximations  of  actual  rotor  behavior.  In  some  instances,  how¬ 
ever,  unsteady  aerodynamic  effects  will  be  an  important  factor.  Therefore, 
the  task  of  including  unsteady  aerodynamic  effects  in  the  Jormal  Mode 
Transient  Analysis  should  receive  further  attention. 
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APPENDIX  I 

AERODYNAMIC  AMD  DYNAMIC  MOMENTS  AND  DERIVATIVES 
IN  THE  RIGID  BLADE  AXIS  SYSTEM 


The  expression  for  the  acceleration  vector  developed  previously  (Eq.  (38)) 
is  employed  herein,  along  with  the  force  and  moment  diagrams  presented  in 
Figures  2  through  1»,  to  determine  expressions  for  the  total  aerodynamic  and 
dynamic  moments  acting  at  a  given  blade  station.  These  moments,  for  con¬ 
venience,  ave  expressed  with  reference  to  the  Rigid  Blade  Axis  System  ("5" 
axis  system) . 


MOMENT  ABOUT 


The  moment  at  a  blade  section  at  r  about  an  axis  which  is  parallel  to 
the  y5  axis  and  which  passes  through  the  elastic  axis  of  the  section 
cam  be  obtained  by  integrating  the  elemental  moment  contributions  due  to 
the  aerodynamic  and  dynamic  forces  acting  on  all  sections  of  the  blade 
outboard  of  the  section  under  consideration.  Thus,  in  Figure  2,  the 
desired  moment  at  station  r  is 


(Ma),5+|MD)>5  .  /  {l mA>ys  —  [(Sa^  +ISld)Z5  +(Sro)l5 
+(sphI25]  +  ISA  ><5  (zseo-Z5„  (r))}df 


/.,T/J%[xs'xsJr)]  - 


a»5(_Z5 -Z5„ 


^rOjjdmdf 


where  (S*)^  and  (SA)Xs  are  the  aerodynamic  forces  per  unit  span  in  the 
Z5  and  *5  directions,  (acting  at  the  elastic  axis)  respectively,  and 
( mA)y  is  the  y5  component  of  the  aerodynamic  moment  per  unit  span. 

The  terries  (SF0)Z  ,  (SL0)Z  ,  (SPR)Z  represent  the  shear  forces  per 
unit  span  in  the5  Zg  direction  introduced  by  a  flap  damper,  a  lag  damper, 
and  the  control  system  pushrod.  The  precise  nature  of  these  forces  is 
described  on  p.  1*9  •  All  quantities  on  the  right-hand  side  of  Eq.  (97)  are 
functions  of  the  dummy  spanwise  variable,  £  ,  unless  explicitly  indicated 

to  be  a  function  of  r  (i.e.,  the  radial  location  of  the  blade  section 
for  which  the  moment  is  being  determined).  Also,  the  differential,  dm  , 
denotes  the  local  blade  mass  per  unit  area  while  the  integral  subscripted  C 
symbolically  indicates  integration  over  the  chord  of  the  blade  section. 


Substituting  Eqs.  (35)  through  (38)  into  Eq.  (97)  and  performing  the  inte¬ 
grations  over  the  chord,  yields  the  following: 


My5'  l  {  (mAly5"[(sA)25  +  tSL0)z5  +  (SPR)Z5+  tSF0)Z5](x5ea-x5#ci(r)) 

+  (S  A  )x5  (zVZ  ^0(r) +.Tm(£  ~  r  ]  {e^  ^e+y10c/e)cose 
-  2^e+y,ocg$e  )^sin  9  -(we+  yIOcg0e) (0sin0+02cos0)  +  v~esine 
+  2ve$cos0  +  (ve+  y  l0c^ )  (0cos  9 -  9%\n9)  +2/3X2  [vecos0 
”  ^ve+  yiocg)®sin  0  -  wesin0-  we  0cosej  +  r  [/sT+/3  (X#+  2 SX2)] 

+  gz5  }d£-jfTn[we  cos  9  +  vesin  9-  we  (r)cos  9  (r) 


-ve(r)sin0(r)]  {-eiV+afty^flsin  9  -  r {X2.2+  2X2S)+  gX5}  d£ 

-  2 {if,  m  ky |o2  0cos2  9  d{  +  J^  my^cos  9  r  ft2d  £  -^my^sinef 
-2ft  [ve  cos  0-^0  sin  0-wesin0-we  ©cosoj-rlftf+a&sj+gxjdf 
r2U^  mkZlo2^sin20dC  (98) 


eft2 


where  the  following  definitions  apply: 


(mass  per  unit  length  in 


direction) 


(i.e.,  blade  section  c.g.  lies  on  the 
y(0  axis  —  this  is  true  for  symmetric 
sections,  and  errors  introduced  for  non- 
symmetric  sections  are  believed  to  be 
negligible) 


(first  moment  of  section  mass  about 
elastic  axis) 


jf  y(0  z(0  dm  =0 

l  yl02  dm  =  mkzl02 


/  t  t 

U.e.,  the  y|Q  axis  is  an  axis  of 
symmetry ) 

(second  moment  of  section  mass  about 
Z;o  axis) 
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/  zJ  dm  s  m  k  2 
Jc  0  yio 


(second  moment  of  section  mass  about 

y»0  8x15 } 


M, 


Equation  (98),  after  nondimensionalization,  becomes 

_  #  ^  J  __  _ 

'y5  *  Jr  {(mA>y5-Sz5l£  "  7')  +  Sx5(25eo-z5Jr))  +  y|0cg(sine 

+  8e  cos 0)  cs0)  x5ly10^=o"2Srn  ( ky|Ocos20+kZ|Osin0)}d£ 

(99) 


where  the  following  definitions  have  been  made: 


Sz5~  ISa,z5+(sfd)25  +  ‘SD  >z 5+<SLD)z  +(SPR); 


(100) 


s>siS*.stsS 


(101) 


(Sn)  *  m 


{  e  +  2  (  v, 


ecos0-ve  0sin0-wesin0  -  we  0cos 6) 


7  ( I +  2S)  -  gx  -  2yi0  0sin0} 

o  eg  j 

(102) 


(  _  ^  X  XX 

l  e/3  +  (  we  +  y,0c  0e )  cos  8  -  2  ( we  +y ^8  e )  8  sin  8 


(S 


o'z,  :  "01 


_  _  XX  x  XX  XX 

-(we  +  Yioc  0eH0sin0  +8cos8)  +7esin0  +  2ve0cos0 


+we  +  y,OJ(0cos0-0sin0)  +  2  /3  cos  ©  -  ( ve+7rocg)  ^  s .  n  0 
"wesin0-we0cos0]  +?  [ /3+/3  ( i  +  28}]  +  qZb) 

(103) 

z5Jr)  =  wec°s  0  +vesin0-we(r)cos  0(r)  -  ve(r)sin5(r) 

(101*) 


‘■'eo  ^ea 
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When  (S0  )x  is  defined  as  in  Eq.  (102),  it  is  recognized  that  upon  sub¬ 
stitution  of  Eq.  (102)  into  Eq.  (99)  certain  products  will  be  neglected 
as  higher  order  in  accordance  with  assumption  10. 

MOMENTS  ABOUT  zs  AND  xs  AXES 

If  a  procedure  is  followed  that  is  exactly  analogous  to  that  followed  in 
developing  the  y 5  moment  and  if  Figures  3  and  ^  are  utilized,  the 
following  expressions  for  the  moments  about  section  axes  parallel  to  the 
Z5  and  x5  axes  can  be  obtained: 

<v  (^A  >z5-hm0)Z5  =  i  { (mA  )Zg+sy5(?  -7 )  -sXg(y5e  -y5 Jr)) 

-  y,0cg  (cos0-0esin  6)  « S0)X5  ^o+0msin20(Kz^-ky|o2)}d^ 

(105) 

MXg  =  (MA)X5  +  (MD\5  =  j[  T{(mA)X5  +SZ5(y5eo-y5eo(r))-syg(z5eo-75M(f)) 

+  y|0cg  [(§D^.yIOc9*o  (cos0-0esin0M§D)y#|y(^,o  (sin0  +  0ecos0)] 

+  m  ( kz2,0  +  k  2o)(-0 4e  -$)+  m  (k2|0  -  kz10)[?  0+2&)sin  20 
■  +  (Q>-/3)cos20]- 2m  [kz2|Ovesin0+ kyi0wecose]  }  d$  (lo6) 


where 


sy5-  (SA)y  +(SD)y  f  (Suo) y.  +  (SFD) y.  +  (SpR)v 


(107) 


(So)y5  =  ~ m  (eS  +  vecos0-  2ve0sin0-(ve  +  y,Ocg)($sin0+0cos0) 

(we+  yioCg^^s'n®“2^e+  yiocg^e^cos^_(we+yiOcg®e)(0cos0-0sin0) 
-2(3  [we  cos  0-  we0sin0  +  ve  sin0  +  (ve  +  yiocg)0cos0]  +  r  ( 8  -  2&k) 
-2yioc.jVe  ~  2/3  [wecos0+(ve  +  y1Ocg)sin0]  -  (1+28)  [(ve+yiocg)cos0 
-wesin0]  +  yiocg0esin0'>gy5}  (108) 


y5co-yt>eo(F)=  vecos0-wesin0-  ve  (f)cos0(r)+-we(r)sin0(r) 


(109) 


s 

I 


Figure  3.  Section  Forces  and  Moments  -  Viev 
Normal  to  Plane. 


39 


Figure  4.  Section  Forces  and  Moments  -  View- 
Normal  to  ytj  -  Zcj  Plane. 


In  Eq.  (107),  ^lo  )y  >  ^fo  V  »  and  ^PR^g  represent  Yg  components 
of  the  lag  damper,  flip  damper,5 and  pushrod  forces,  respectively.  These 
are  given  in  Appendix  III. 

SPANWISE  DERIVATIVES  OF  MOMENTS 


The  following  expressions,  which  represent  the  first  and  second  derivatives 
with  respect  to  7  of  Eqs.  (99)  and  (105)  and  the  corresponding  first 
derivative  of  Eq.  (1C6),  are  also  required  in  the  derivation  of  the  final 
modal  equations  of  motion: 

"v  r  ,f)5*sK 

2  2  \ 

+20  fn  ( ky|Ocos20  +  kZ|Osin20) 

(lio) 

%5=  “SZ5~  z'l0  f/(Sols  dC  +  Z5ea  (s0)X5- [y.Ocg  Sir,0(§D)x5,y|Oc^O  +  (™4, 

-20m  (ky|Ocos20+kz^osin2e)] 
(ill) 

{“ §y5+ ho lF)  sx5}df +y,ocgcos 0  (§0^,y10cij=  0  "  (rH5 


x  /-2  -2  \ 

-m0sin20  \kzl0“  kyio/ 
(112) 


Mz5=Sy5-y5eaSx5+y5eo  f?  sx5dC  +  [yiocgcos®(SD)x9  y(0  :0-(mA)Z5 

CQ 

-m0sin  20(kz,o-kyio)] 

(113) 

“v  J-  §2sd^  +  *L  Jf  \d*  " M.9  -  y.ocg  [(5D)l9  cos® 

-(SD)y5S«n©]yi0cg=o  +m(kZ|0+  kylo)(0+4X+/S )-  m(ky|0 

~  * zio )  [  2  ( 1  +2  S)sin  2  9  +  (0e  - P)cos  20J+  2  m  kz |Q ve  sin  0  +  2 m  ky^ecos0 

(114) 


where  the  angle  ©  is  given  by  Eq.  (45). 
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APPENDIX  II 

AERODYNAMIC  SHEAR  FORCES  AND  MOMENTS 


To  determine  the  aerodynamic  shear  forces  and  moments  acting  on  the  rotor 
blade,  a  hlade  element  approach  was  followed  that  utilized  the  quasi-steady 
aerodynamic  theory  as  developed  in  Reference  5*  Quasi-steady  theory  neg¬ 
lects  the  effect  of  the  wake  generated  by  the  blade,  and  its  use  is  dic¬ 
tated  by  the  lack  of  a  generalized,  variable-inflow  theory  applicable  to 
rotating  wings  in  forward  flight.  In  the  analysis  herein,  the  apparent- 
mass  terms  appearing  in  quasi-steady  theory  have  also  been  neglected  as 
second  order.  The  primary  differences  between  the  resulting  quasi-steady 
aerodynamic  expressions  and  those  predicted  by  classical  steady  aerodynamic 
theory  are  that  (l)  the  aerodynamic  forces  are  determined  by  the  velocity 
components  at  the  75#  chord  point  and  (2)  an  aerodynamic  damping  moment  in 
pitch,  fhj  ,  is  predicted.  The  use  of  the  blade  element  approach  implies 
that  radial  flow  effects  are  also  neglected. 

When  the  aerodynamic  shear  forces  and  moments  are  derived,  reference  will 
be  made  to  Figure  5, which  shows  a  cross  section  of  the  blade  as  well  as 
the  relative  air  velocity  components  UT  and  Up  .  The  tangential 
velocity  vector  UT  is  parallel  to  the  plane  of  rotation  of  the  rotor 
and  normal  to  the  local  X,0  axis  while  Up  is  normal  to  both  the 
local  X|0  axis  and  UT  .  UT  and  Up  thus  lie  in  the  plane  formed 
by  the  blade  section  axes,  y[0  and  zio  •  Within  the  assumptions  of 
small  elastic  rotation  of  the  blade  section,  the  incidence  angle  €  is 
equal  to  the  local  total  pitch  angle  ©  .  The  aerodynamic  force  and 

mcment  components  in  the  "lO"  axis  system  are  then 


<SA>zI0  = 

(SAl  = 

amo 


(  /  cos<£  +  dsin<£)cos© 


(^sin<£  ~dcos<£)cos® 


(SJ 


A'x 


10 


10 


-  (/sin <j)  -  dcos£)sin© 

(115) 

+  ( /  co$(p  +  dsin<£)sin© 

(116) 

0  '  (117) 

0  (118) 

0  (119) 


(mA)X[o=  mc/4  +  Ucostp  +  dsin<£)y10c/^cos©-  (£sin<£  -  dcos<£  J^o^sinQ+n^ 


1»2 


(120) 


(126) 


In  the  above,  J2.  is  the  local  lift,  d  the  local,  drag,  mc/4  the 
local  pitching  moment  about  the  25$  chord  predicted  by  steady-state  aero¬ 
dynamic  theory,  and  the  damping  moment  predicted  by  quasi-steady 

theory.  With  use  of  conventional  aerodynamic  notation,  the  following  ex¬ 
pressions  are  readily  derived: 


(130) 


In  Eq.  (130)  the  quantity  a0  represents  the  nondimens ional  distance 
(in  terms  of  50$  chord  lengths)  from  che  mid-chord  to  the  pitching  axis  of 
the  blade  section.  The  quantity  a0  is  positive  when  the  pitch  axis  is 
downstream  of  the  mid-chord.  For  a  rotor  blade,  the  downstream  direction 
depends  on  whether  the  ^lade  is  in  conventional  or  reverse  flow.  Hence,  o0 
is 


a0 


2_ 

c 


2 


(conventional  flow) 


1 

2 


(reverse  flow) 


(131) 


^5 


The  aerodynamic  coefficients  in  Eqs.  (127)  to  (130)  are  functions  of  the 
section  angle  of  attack  and  Mach  number,  as  determined  by  the  velocity 
components  at  the  75#  chord  point.  From  Figure  5>  the  angle  of  attack  is 

Or  ”  ®  +  Ian"'  77  (132) 

Ut 


The  relative  air  velocity  components  are  given  by 

777  =  ( V^io).3. c cos ®  -  (Vzio)^r  sin®  (133) 

SIR  4  4  c 


UP 

SIR 


cos® 


K),c  sin® 


(13»0 


Or.,,  expressing  the  "10"  axis  velocity  components  in  terms  of  "5"  axis 
components,  Eqs.  (133)  and  (131*)  are 

-Hi  =  (Tys)3c  +  (%),  (XlSin@  "  X2C0S  ®) 

ilR  4l  4 

(135) 

^  =  -  Nic  +  (  X|C°S@  +  X2Si"®) 

(136) 


By  using  Eq.  (29)  to  determine  ^z^z/4  c  »  etc.,  in  terms  of  "1"  axis 
velocity  components ,  by  neglecting  higher  order  terms, and  by  noting  that 


can  be  written  more  familiarly  as 

V0x  =  ,  vcoso?  = 

°XI  ilR 


(137) 


Eqs.  (135)  and  (136)  can  be  rewritten  as 


=  ^  f(l  “  -|2)  s?n  »//  +  Scos^/J  +  1  +  ?  (i- ^  +  |)  +  7ecos B 

+  y,o^c)|sine  -  (4  +  7,„?cle)  sine  -  (  we+  y,0}9«)|cose 
'  V6'  “  ^fWjCose  +  (7e  +7i0|c)sin  9j 

+{ WgS in e  -  Vcos q )  [-/.(cos  *  -  8  sin  *)  -  X./}-  yl0  Cos9  J 
-0'|vasin0  +  wecos  A1  cos  ^ 


(139) 


ftR  Xs  ('"4")  ~ /i/?(C0S^  ~  $sin^)  -  (we+ yi0_  JeJCos0 

+  ^we+  yio^c0ej^sin0  -  vesin0  -  ^ve+  ytOicj0cose 

-/3pVe  + y|Oicjcos0  -  We  sinfi  j  -  r/3 

+(we'cose  +  Vg'sinfi)  J^-r.  (cos*-  Ssin*)  -  xsfl-  cosfll 
-^COSI//  (veCOS0-w  cinn^  C  J 


(1^0) 


Now,  inasmuch  as  the  aerodynamic  terms  (^a)x  (m*)  y«  (rrw) , 

Thus^the^inil  CO°l'Date  faCt°rS  ’  th6y  Can  be  °eSlected  as  higher  order. 

us,  the  final  aerodynamic  force  and  moment  expressions  are  (after  sub 
stitution  of  Eqs.  (127)  to  (130)  into  (121)  to  (125)) 

^*5  =  i[w)z  ($gf(cxCos$  +  cdsin*)  (m) 


^)'s  ‘  -  °d«>s*) 


(ll»2) 


+7,0^  ($f  [oose  (V05*  +  C,iSi"<f'' 

-sin©  ^c(sin<#p  -  cdcos^)]j  ( 


(H>6) 


The  following  pitch  ang^e  relationships  are  also  required  in  the  aerody 
nanic  expressions : 


Total  section  pitch  single 

®  =  5  +  0e  *  e0  +  eB  +  et 

Built-in  pitch  angle 

0B  =  6,  (5  +  7-0.75) 


.  (U»7) 


(lh8) 


Control  system  pitch  angle 


a0  =  e,  -  A,  cost  -  B,  Sint  -  ton  83  [£ +(wA] 


APPENDIX  III 

FLAP  DAMPER,  LAG  DAMPER,  AND  PUSHROD  FORCES 
FLAP  DAMPER  FORCE 

in  this  analysis ,  it  is  assumed  that  the  flap  damper  force  is  in  the  Z4 
direction  and  is  located  at  a  distance  rpo  from  the  flapping  hinge. 
The  damper  thus  provides  a  pure  moment  about  the  flapping  axis  y4 
which  can  be  expressed  as 


Mp0=  rFD  cos/3  (spd)z41t 


(150) 


where  ( SF0)T 
damper 


is  the  total  damper  force  in  pounds.  For  a  pure  viscous 


(sFD)z4,r  cfd(t7£asl)(;| 


(151) 


c  * 

where  ''fd  is  the  line*  ”  damping  coefficient  of  the  damper  in  units 
of  lb-sec/ft  and  the  bracketed  expression  represents  the  approximate 
velocity  of  the  blade  in  the  Z4  direction  at  the  point  of  attachment 
of  the  lag  damper.  By  nondimens ionali zing ,  and  if  Eqs .  (12),  (36)  and  (37) 
are  used,  and  if  higher  order  terms  are  neglected,  Eq.  (15I)  becomes 


<5ro):4,T  =  --^ri^  +  "« cose 


—  we  0sin0  +ve  sin0  + 


ve  0  cosS  ] 


(152) 


Resolving  the  above  force  into  "5"  axis  components  and  neglecting  higher 
order  terms  yields  the  following: 


(Sfo)z5;t  ~(sfd)z 4,t 


(153) 


(SFD)y  T«0 

*R  '  ' 


(1510 


i+9 


To  obtain  the  flap  damper  force  per  unit  span  >  it  is  -ssumed 
that  the  total  damper  force  is  applied  uniformly  over  the  infinitesimal 
interval  A  .  Thus , 


^  .or(7FD-A)srs(TFD+A) 
0  for  (Tro+4)<T  <(7FD-|-) 


(155) 


LAG.  DAMPER  FORCE 

The  required  lag,  damper  force  per  unit  span,  (^Lo)y^  ,  can  be  obtained 

from  an  analysis  similar  to  that  for  the  flap  damper  assuming,  however, 

£. that,  the  damper  force  acts  in  the  yr  direction.  The  following  result  is 
obtained: 


A 

^  rLD 

A 
”  2 

y  ^  (rLD+- 

*■) 

(SL°)y5  =  < 

0  for  |VL0  +• 

A  \ 

L  A\ 

t) 

<7<( 

(156) 


(157) 


where 


(Sld)V= 


(5  LD)ys,T  a  CtD* 
m0ftzR2  m0iiR 


_  *  JL  _  x 

78  +  ve  cos6- v e  0sin0- 

_x  _  XT 

wc  sin©  -we  0cos0 

J  LD 


(158) 
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PUSHROD  FORCE 


Rotor  blades  are  generally  restrained  in  pitch  (i.e.,  the  motion  a?;out  the 
*5  axis)  by  the  moment  exerted  about  the  xg  axis  due  to  the  pushrod 
o f  the  control  system.  For  conventional  rotors,  the  pushrod  force  passes 
through  the  y5  axis  and  is  parallel  to  the  Z4  axis.  A  pure  pitching 
moment  is  thus  provided.  For  rotors  having  pitch- flap  coupling,  the  push- 
rod  is  displaced  from  the  y5  axis  in  the  direction  of  XR  by  a  dis¬ 
tance  rPR  •  In  this  position, the  force  of  the  pushrod  exerts  a  moment 
about  the  y5  .axis  and  thus  influences  the  flapping  and  bending  motions 
of  the  blade  as  well  as  the  torsional  motion. 


Pushrod  effects  in  the  torsional  equation  are  introduced  by  representing 
the  torsional  deformation  by  a  series  of  uncoupled  modes  determined  on  the 
basis  of  appropriate  roo*.  boundary  conditions  (i.e.,  fixed  or  partially 
fixed).  To  determine  the  effect  of  the  pushrod  force  in  the  blade  flap¬ 
ping  and  bending  equations,  an  expression  for  the  pushrod  force  is  required. 
With  the  representation  of  the  torsional  response  of  the  blade  in  terms  of 
uncoupled  modes,  the  required  force  is  approximately  given  by 


(spr) 


[g->  fys/ ggiL, 


(159) 


(sm  )w  -  0 


(160) 


where  ■  X 5PR  is  the  distance  from  the  *5  axis  to  the  pushrod.  The 

corresponding  pushrod  force  per  unit  span  is 


(s>R) 


'or(7PR-A 


)S?S  ( 


rPR+if 


0  for  ^  rpR+  2  )  <  r<(rPR~  2  ) 


(l6l) 
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APPENDIX  IV 
MODAL  CONSTANTS 


various  constants  appearing  in  Eqs.  (80),  (86),  (89)  >  (9l)  and  (98) 
given  "below: 


c 

(162) 

0 

m0fl0Ri 

= 

/0rTffi*>idr 

(163) 

c*i,p 

= 

(16U) 

s 

fTr-  l 

Jo.  mV«i'T<lf 

(165) 

C5. 

Di>P 

= 

/0s>»iVd? 

(166) 

—> 

= 

(167) 

C«i.P 

= 

C*J?  '*  dF 

(168) 

CV 

= 

e'/o?-,[ -  f 

(169) 

r> 

ffT 

HOi 

J0  myw(dr 

(170) 

CI2| 

s 

X  "%j7d? 

(171) 

C1V 

/•rT 

Ibj 

9.  mr  v  (7-0.75)dr 
>*'o  "i 

(172) 

C|4p 

= 

X  ^pd? 

(173) 

(•ft 

CI5P 

z 

Jo  »*,* 

(17U) 

C,8P,p' 

= 

X  "’yyrdF 

(175) 
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C2°p,p'  = 

(176) 

C2'p,p' = 

el\  ^,TS/v(f-7]  dj  dr 

(177) 

C22i,p  = 

A't  rfr 

JoVr  "%id/;dr 

(178) 

C23p  = 

-  2V*  “*A^]d?  +  [as/£/?  Xd^ 

-  2  A5A^/ncw  ~  A?Acwm^]dr 

(179) 

s  ■ 

j[  Vvpte^WdC^e/mT-e^mV+m)]^  +  jf  >vp[A«A^ 

°ew - 

fficwCdC 

-2 A e/  mcwr  - AeA  (m^r  +  fnCT()]dr 

(180) 

S?6ltP  = 

aa1J^p(EXlX,)"dr  S-28,^[(Ei,^'  + 

(181) 

C27lfP  = 

S'0rJ&iriJ *  tEI»R  +  n»4ldf 

(182) 

C33p  = 

®'[  /0r'ffiv2<17  -  a75C«J 

(183) 

C34pip/  = 

eJ0  ffiVy(f" 0,751  dr 

(184) 

C36j  = 

jTVkf^^dr 

(185) 

c37j 

^XV5*  /f  T  s'? d?]  “f 

(186) 

C38i/,j  = 

2X  V»iy»“(EIl'EI>ldf 

(187) 

C39.  .  - 

33J.P,p' 

‘  2  X  >&/- Pr-p' (fI'  “n»)d7 

(188) 

C44.  = 

X 

(189) 
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(205) 


>8,'  l 


rT 

X. '  =  4,  %  0.75)  dr 

69i,i'  h  m  'w|  7W,  y,0cs  dr 


o 

-'J 

o 

■o 

H 

r~ 

J  m  y  y  dr 
*o  %  yiocg 

V72p  = 

, 

J  rv  my,0  dr 

0  vp  °cg 

rh 

C73p  = 

J  mYv  y,o  7dr 

o  VP  ucg 

r?T 

^74  .  = 
P,J  ' 

/t  /  /*fT 

n 

-Ni 
JJ 1 

5» 

n 

J0  yvPV  4  mn,ld?<r 

/- 

'C76p  =  * 

C7W=, 

/T  /•*T 

4  V  4  m  v*1* 

/T_ 

c  = 

79p,p'  ' 

/  m  y  y  '  y  dr 

0  %  ^  y,0c, 

m  Ye .  kz  V  -  0.75)dr 


ffT_ 

Cso j  -^i  J0  m 

/T_ 

XHmV«Wdr 


-  f'T- 

=  J  m. 

i  ■'r\ 


77  2 


S?  .=  J  m'  k,  "  yft  y  7  dr 
■  82P,j  •'o  zio  rB\  \  u' 


(206) 

(207) 

(208) 

(209) 

(210) 

(211) 

(212) 

(213) 

(214.) 

(215) 

(216) 

(217) 

(218) 
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j T 

c„  zJ  myl0  yJ dF 

83i,J  'o  ®J  l0cg  wi 

(219) 

/ 

c«u ■  4  m  V*,  Vdr 

(220) 

c»p,i  = 

(221) 

c.6iil  ■  r\<krw-? + AsAj0rc>(«+^]dr 

(222) 

c%,i  ■  /0Vv:[-54'Wdf+A\i'  m»M]dF 

°cw 

(223) 

c«1.i  *  4\<[-5*/r,’iSfdf+A\.-(ircJ,ic»fdf]dr 

« 

(22U) 

cTj 

c»op  ■  4  *  v-;d7 

(225) 

-C\(mK)6T 

(226) 

/t 

c92j  =J0  ^V7 

(227) 

/* _ 

c93j  =  J0  myiocgyejrdr" 

(228) 

Jr 

c9^--  «i  J0  m?ioct>'8)7(r-a75)dr  ' 

(229) 

CM.  =  ®!  4  'mV9|  yw|(r"°'75)dF 

(230) 

c,e  *  !j\f 

(231) 

C97  =  l  Tfi-ry,ocf 

(232) 
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(233) 


C"l  D  = 

•*  P 

l  **!*,>»!* 

(23U) 

C'00p,i = 

/o'my»Py8lfdF 

(235) 

C"°i,i  = 

/  m  k_2  x  y  '  dr 
:o  zio^j  % 

(236) 

C|lli,P,j  " 

A 

i  m  V«i  <  yvB  df 

(237) 

•CmKy«i  r"/w|/dr 

(238) 

A 

J0  Vdr 

(239) 

c"<p,i  = 

yvpff-°-75)dr 
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